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I Abstract 

. The instanton solution for the forced Burgers equation is found. This solution de- 

, scribes the exponential tail of the probability distribution function of velocity differences 

I in the region where shock waves are absent. The results agree with the one found re- 

cently by Polyakov, who used the operator product conjecture. If this conjecture is 
true, then our WKB asymptotics of the Wyld functional integral is exact to all orders 
l/^ I of perturbation expansion around the instanton solution. We explicitly checked this in 

■ the first order. We also generalized our solution for the arbitrary dimension of Burg- 

ers (=KPZ) equation. As a result we found the angular dependence of the velocity 

difference PDF. 
Oh! 
(D . 

^ . There are two complementary views of the turbulence problem. One could regard it 

as kinetics, in which case the time dependence of velocity probability distribution function 
' (PDF) must be studied. The Wyld functional integral describes the correlations functions of 
5^ ; the velocity field in this picture. 

Another view is the Hopf (or Fokker-Planck) approach, where the equal time PDF is 
studied. For the random force distributed as white noise in time, the closed functional equa- 
tions (Fokker-Planck equation) can be derived. In the case of thermal noise the Boltzmann 
distribution can be derived as an asymptotic solution of this equation. 

One of the authors 0] reduced the Hopf equations for the full Navier-Stokes equation to 
the one dimensional functional equation (loop equation). The WKB solutions of this equation 
were studied, leading to the area law for the velocity circulation PDF. 

In the recent paper by Polyakov |^ a similar method of solving the randomly driven 
Burgers equation was proposed. It reduced the problem of computations of their correlation 
functions to the solution of a certain partial differential equation. This equation for the 
velocity difference PDF can be explicitly solved. 

The derivation of the Polyakov equation was based on the conjecture of the existence of 
the operator product expansion. 
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On the other hand, it was recently understood in the paper by Falkovich, Kolokolov, 
Lebedev and one of us that the Wyld functional integral provides the general method of 
computation of the PDF tails. One has to find an instanton - the minimum of the action in 
the Wyld functional integral. 

This instanton is not the same as the solution of the classical equation in the usual sense. 
The force is present, and it acts in a self-consistent way, as required by minimization of the 
action. This force is no longer random, it is adjusted to provide the large fluctuation of 
velocity field under consideration. 

It was shown in that the PDF for the passive scalar advection in the Gaussian velocity 
field is asymptotically described by an instanton, with the spatially homogeneous strain. 

Here we find instanton in the Burgers equation, in the presence of finite viscosity. The 
result which we obtain in the turbulent limit (vanishing viscosity) coincides with that of 
Polyakov which gives an indirect confirmation to his OPE conjecture. 

We start with the randomly driven Burgers equation. 

Ut + uu^ - vu^^ = f{x, t) (1) 
where the force f{x,t) is a gaussian random field with a pair correlation function 

{f{x,t)f{y,t')) = 6{t-t')K{x-y) (2) 
The Wyld functional integral has the following form (see, for example, |^). 

I Vfexp (^-^Jdxdydt f{x,t)D{x-y)f{y,t)^ = (3) 

j VfVu 5\ut + uu^ - uu^^ - f] exp (^-^ J dxdydt f{x, t)D{x - y)f{y, t)^ = 

J VfVuVnexp (^i J dxdt^{ut + uu^ - vu^^ ~ dxdydt f{x, t)D{x - y)f{y, t)^ = 

J VuVfi exp J dxdtfi{ut + uUx — vUxx) ~ 2 J ^^^y^^ /i(x,t)/t(a; — y)fi{y,t)^ 

Here we started with the obvious functional integral for the force, where D is the function 
inverse to k. To change the variable of integration from / to m we inserted into the functional 
integral the identity 

M = j Vu 5[ut + uu^ - vu^x - /] (4) 

where M is just a number. Dropping that number as an unimportant constant, we removed 
the delta function at the expense of introducing a "conjugated" variable n and evaluated the 
integral over / to arrive at the final expression in (^. 

It is not so obvious that the integral in (^ is equal to a pure number because of the 
determinant det arising in its computation. It is possible, however, to prove (^ using 
the causality argument (see 

Thus the initial problem of computing the correlation functions of Burgers equation is 
reduced to the field theory with the action 

S = -I J fi{ut + UUx ~uux:r) + ^ J dxdydt fi{x,t)K{x -y)fi{y,t) (5) 
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Here we are going to study the correlation function 



(exp{Ao(u(po/2) -u(-po/2))}) = J VuVfiexp {Xo{u{po/2) - u{-po/2)) - S} (6) 

whose Laplace Transform gives us the two point probability distribution (see 0). We will 
often refer to the expression we have in the exponential as the action 

Sx, = S- XoHpo/2) - u{-po/2)) (7) 

There are no general methods to compute the functional integrals like (^ exactly. The most 
straightforward approach would be to expand the exponential in functional integral in powers 
of the nonlinear term puUx- By doing so we will just reproduce the well known Wyld's diagram 
technique (see 0). The attempts to use this technique to describe turbulence always failed 
because we are interested in the limit i/ ^ when nonlinear term dominates the functional 
integral. The absence of a large parameter makes the task of computing this functional 
integral by using perturbation theory hopeless. 

Nevertheless, if we are interested in computing the large Aq behavior of the correlation 
function in (^, we can use Aq itself as a large parameter. Then the integral will be dominated 
by its saddle point, or by the solutions of the equations of motion for the action (|^). All 
we have to do is to find those solutions and compute the value of the action Sinst on those 
solutions. The answer will be given by 

(exp {Ao(^/(po/2) - «(-po/2))} ) = exp {-S,^,, (Aq) + S,^,, (0)} . (8) 

If we want, we can then further expand the integral in powers of I/Aq using the perturbation 
theory around those solutions. We will call this method WKB approximation and the solutions 
instantons using the names borrowed from quantum field theory. 

To that effect, let us write down the equations of motion corresponding to the action (|^. 
They are 

Ut + uu^ - vu^x = -I I dy k{x - y)p{y) (9) 



pt + up, + upxx = -I Ao - f ) - + f ) } '^(^) (10) 



To solve these equations, let us first notice that the only role the right hand side of ([10|) 
plays is giving the field p a finite discontinuity at t = 0. It is also easy to see that p{f) = for 
t > 0. This is because p feels a negative viscosity so any solution which is nonzero at t > 
will become singulaiQ. Thus the field p can be evaluated at t = — to be 

^{t = -^)=%\Ab[x-^-^~b{x^^-^ (11) 
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while it is zero at all later moments of time. It is therefore convenient to speak of the field p 
propagating backwards in time starting from its initial value given by (|TT|). 

^ Alternatively, one can argue that the integrals in (^) are defined only for i < 0. Those arguments use 
a striking similarity between and a Feynman path integral for a quantum mechanical system with the 
coordinates u and momenta /i to define (g) as a wave function in the momentum representation. Then the 



conditions (11) become obvious. 
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If we try to propagate ( [Tl| ) back in time, we discover that we have to deal with two 
phenomena governed by the second and third terms of ([l0|) . One of them is just a motion of 
the initial conditions as dictated by the velocity u. The other is the "smearing" of the initial 
delta function distributions in ( pT]) due to the viscosity. 

However, it can be shown by a direct computation that the smearing does not change 
the value of the action on the instanton as long as the viscosity is not very large. We will 
construct a solution which takes into account the viscosity in the end of the paper. For now 
we will just drop the viscosity term to arrive at a simplified equation 



fit + M/ij 







(12) 



Since all this equation can do is moving the 6 function-like singularities around (and changing 
their heights by compressing them), it is clear that the solution of (|l^) with the boundary 
conditions given by (|TT|) is just 



fi{t) =iX{t){6 [x 



Pity 



X 



Pity 



(13) 



with the boundary conditions 



A(0) = Ao, p(0)=po 



(14) 



Now let us leave the equation (|TOD for a while and study (^. A natural thing to do is to 
substitute ([l3|) into the right hand side of (|^). 
We obtain 



Ut + uux — jyu^ 



X < K { X 



- — K [ X -\ 

To proceed further we need to know n. Let us assume, following 
varying even function of x which behaves as 



k{x) ^ k{0) 




= L 



(15) 

that k{x) is a slow 
(16) 



and quickly turns into zero when |x| ^ L. The interval L characterizes the range of the 
random force and we will work only there, that is we suppose that p also lies within this 
interval. It is clear then that the contribution to the action (^ comes only from the interval 
L (compare with (0)). So we do not have to know the velocity beyond that interval. There 
we use (|T6|) to obtain 

Ut + UUx - Z/M^^ = XpKQX (17) 

Notice that k(0) dropped out. 

(p!?!) is a Burgers equation with a linear force. It is easy to solve such an equation. We 
have to look for the solution in terms of a linear function 



which leads to 



u{x, t) = a(t)x 
da 



dt 



a 



kqXp 



(18) 
(19) 
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Notice that the viscosity term did not contribute. That does not mean that the viscosity is 
not important at alL For a; 3> p the force in the ([l5|) becomes zero and the viscosity there is 
important to make the velocity go to zero at the infinity. However, in the region x o: p which 
is the one we study the viscosity term can be dropped. 

Now we can use (0) and (|TBp to solve (|T^. A direct substitution leads to 



H 
dp 

11= 



Xa (20) 



These can be solved in terms of the function 

R{t) = exp (^J^ dt'a{t')^ 

to give 



A = XoR (21) 
p = poR 



while R itself satisfies, by virtue of ([T9|), the equation 

d^R 



dt^ 



KopoXoR' (22) 



The last equation has to be solved with the boundary condition i?(— oo) = otherwise the 
action ([7p will not be finite. The solution is given by 

R = ^ (23) 

1 - 

So we have found the instanton solution for the equations (|^), (|10|). Notice that it is the only 
solution of the equations of motion, so we do not have to sum over different instantons. 

Now it is a matter of a simple computation to find the action on the instanton. We collect 
everything together and substitute (^3]), (pi]), (0), and ([18|) back to (^ to get 



5inst = -^(AoPo)^ (24) 

while the correlation function we were studying is 

( exp {Ao (n(po/2) - u{-po/2))} ) = exp (AoPo)^) (25) 

This is the same answer as the one obtained in 0. We want to emphasize, however, that we 
obtained it without any conjectures and only as an asymptotics for |Ao | ^ 1. 

We are not going to discuss the physical implications of (p5| ) referring instead to the papers 
i and i. 
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Now we return to the question of why we can drop the viscosity in (|T0|). Namely, we 
just construct the solution of ([l0|) with the viscosity. To do that, it is convenient to Fourier 
transform it (taking into account that the velocity m is a linear function of x). 

^ " (P/^l/')) - ^ = (26) 



Then the solution of (|26|) can be found as a direct generalization of (|13D. 



/i(p) oc \{t) sin f 1 exp (-/9(t)p') (27) 



Here we had to introduce the new variable (3{t) which measures the speed of smearing of the 
solution with the evident initial condition /3(0) = 0. Substituting (^) into (|^) we reproduce 
the equations (^TD for p and A with the additional equation for /? 

pt-2aP + iy = (28) 

Pit) = (29) 



with the solution 



where uj = sj {kqPqXq) /2. 

The smearing, however, has no influence whatsoever on the velocity. To see that, we 
substitute (^) into (^ and arrive back at (|l^). In other words, the variable a still satisfies 
the same equation (|T^. A simple argument given below shows that the value of the instanton 
action depends only on the final value of the velocity of the instanton solution which, as we 
just showed, does not depend on the viscosity. 

We must remember, though, that we should not allow /i to spread beyond L interval, or 
more precisely L > a/^. (3 can become arbitrarily large for large negative times, but the 
characteristic time interval which contributed to the computation of the instanton action is 

tinst = - (30) 
LJ 

So all we have to do is to make sure that /^(tinst) < L'^ or 

> - (31) 



This is the condition which viscosity must satisfy for (^5|) to be correct and independent of 
viscosity. 



So far we cannot claim that (^5]) is an exact answer. It is just a leading asymptotic if |Ao| 
is a large number. It might be important to estimate the next order contribution to the ( P3| ) 
especially in view of the claim made in that (p5[) is actually exact. 

To do that it is convenient to introduce the quantity 

^ , / f f I /oxu\ ST^uVp (m(po/2) - m(-po/2)) exp (-^Aq) 

— log expjAo M po/2 -M -po/2 } = — — 32 

dAo j VuVp exp (-Saq) 
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It is easy to expand this quantity around the instanton solution. Writing u = Winst + u, 



A* = /^inst + /i we arrive for (|32D at 



'P0\ ( P0\ , /I'MP/i (m(Po/2) - M(-po/2))exp {-Sao (ttinst + M,/^mst + /i)} 

''^inst 1 TT ) "~ "^inst 77 



2 ) \ 2 J J VuVfi exp {-^Ao (Winst + U, /iinst + p)} 

(33) 

The first term of the expression corresponds to the instanton contribution to the action. 
We see that it depends only on the value of the instanton solution at t = 0, justifying our 
argument that the viscosity does not contribute to the answer. We can easily compute this 
term by substituting the known instanton solution to get 



y^Po (34) 

which is of course compatible with (p5|). 

Let us proceed to estimate the second term in (|33|) . It is not difficult to see that if we 
expand ^Aq (winst + u, Pinst + p) in powers of u and /i up to a second order, the contribution 
from that will be zero as S will be an even function of u and jl. However, due to the presence 
of the third order terms in the action, it is not clear if the higher order terms are also zero. 
If they are nonzero, then a dimensional argument shows that the next term in (|^) is of the 
order (1/Ao)^. Really, then the expanded action looks like 

S{Uinst + U, /iinst + p) = 5'inst + 5'quad + "S'int = 'S'inst ^ ^ J pi^t ~ l^Uxx) + (35) 
^ J dxdydt fj,{x, t)K{x - y)fL{y, - i J (uUinst) - i J PmstUUx + Sint 



'int 



'S'quad meaning the quadratic part of the expansion and S'int being the interaction term, Su 
—I J jluux- We need to take into account the interaction in the first order of perturbation 
theory. That is we need to evaluate 

{U (po) S'int) (36) 

understanding the average in the sense of the action Squad- While the "honest" computation 
would require the computation of the Green's functions for that action, we can estimate the 
value of (|36D by noting that Winst oc \/Ao, /iinst Aq and 



S'quad OC Ao + Y Ao /iM (37) 

Therefore the typical fiuctuations are 5u oc and 5jl (x 1. The expression (^) involves 



three u and one /i, so we could expect it to be of the order I/Aq . However, there is also an 
integration over t involved. That integration gives us l/cu oc I/a/Aq. Altogether we arrive at 
1/A^. 

In other words, it could be expected that the correction to the answer has the form 



>xp {Ao (m(po/2) - m(-po/2))} ) oc exp ( (AoPo)^ + (38) 



^2kq ^ ^ ^1 ^ const \ 
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However, this question requires much deeper investigations which we leave for future work. 
The formula (^) is just a dimensional estimation. The Green's functions for Squad have to 
be constructed for the corrections to the answer p5| ) to be computed reliably. That has not 
been done yet. 

The analysis of this paper can easily be extended for the case of more than one dimension. 
The analog of the equations (|) and ([T0|) will be 

Ou I d \ f 

'dt ^ l^-' ar" ) ~ ^ ~^ J '^'^^^ ~ '^^^^^y^ ^^^^ 

The solution of these equations is a direct generalization of ([13D, or 

pM = ^ m - ^) - ^ + } (41) 

with 

A,,(0) = Ao,,, p^{0) = Poi (42) 

The further progress depends on the tensorial structure of which is just a correlation 
function 

{Mx,t)f,iy,t')) = K,jix - y)Sit - t') (43) 

A natural thing to assume would be that the force is a gradient of something, that is = c?j$, 
in which case 

Kij{x) ^ Kij{0) - y (x^5ij + 2xiXj) (44) 
A direct generalization of the velocity ansatz is 

Ui = aijXj (45) 
and the equations (|19]) and (|20|) turn into 

^ = A,^.. (46) 
dpi 

IF = ""'^P^ 

a can actually be eliminated from those equations to give us an analog of 

d'^X, 



[piX^ + 2XAkPk) (47) 



dt^ 

d"^ Pi I \ 1 o \ 
= 1^0 [XiP + 2piXkPk 
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While a general solution of those equations is rather difficult to find, it is possible to find 
the action on the solution by analyzing the corresponding Hamilton- Jacobi equation. To do 
that, we note that the equations (P7| ) are hamiltonian with the hamilton function 

The (time independent) instanton action 5* clearly satisfies the equation |^ 

dSdS (^2^2 ^ 2(A.p.f ) = (49) 



dXi dpi 2 

By rescaling the time and the variables p and A in (^) we can show that the action S has 
the following initial condition dependence 

^ = /y(PoAo)V(cosvp) (50) 

where ip is the angle between the vectors of initial conditions Aoi and poi- This ansatz coincides 
with the one dimensional answer ( ^If ) up to a nontrivial function of the angle / (cos (p) which 
we would like to determine. Plugging the ansatz into the Hamilton- Jacobi equation we obtain 
the equation for /(cosy?) 

^zf + 3//'(l - z') + f\-z + z') = l + 2z^ (51) 

where z = cos((y9). This first order differential equation has to be solved with the boundary 
condition 

/(I) - ^ (52) 

which follows directly from the equation (^) but also can be computed by solving the equation 
of motion for ip = 0. We can find the function / as a series in powers of 1 — z. It turns out 
there are two solutions 



2 VS + Vn , , 5v/33-61 2^ 
f(z) = —pz 1 — 2; H 7 — — =r{l — z) +... 53 



2 yn-VS , 5^33 + 61 ,^ o 

= —j=^ : (1 - ^) 7 — p 1=\{^ - -2) + • • • 

W3 4 32 f3V3 + 2v^r 



The equation (|51| ) does not tell us which of these two to choose. We have to match (|53| ) 
with the solution of (^Tf). Those equations cannot be solved in general, but there is a way to 
find their solution if the angle ip is close to zero which should be enough to determine /'(I) 
and therefore to choose the right action. 

^We would like to note that the equation (^9|) follows from the master equation of |^ if the viscosity is 
completely neglected. The authors are grateful to A. Polyakov for pointing that out. 
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To do that, we note that the motion represented by ( pT]) is essentially two dimensional, 
with all the motion confined to the Aq, po plane. Then we represent A as a two-vector (Ai, A2), 
while p = (Ai, — A2). The equations (^) turn into (we choose the units where k,q = 1) 

'^^Ai ^3 2 



3Af - AiA^ (54) 



— 3A2 + A2A^ 



We choose the boundary conditions Aq = 1, po = 1. If = then A2 = while Ai satisfies 
the equation 



d'Xi 
IF 

hence 



3Af (55) 



Now if is a small number then A2 ^ Ai and it satisfies the approximate equation 

dt^ {l-utf 

with the solution 



^ {l-uty 6 
In particular, for t = 0, 

d\ogX2 -73 ±v^ 



C -3 ±733 

a = (58) 



dt 2V2 

That last quantity can also be evaluated if we know the action S. For t = we obtain 



(59) 



dpi dS 3Xi f'{z) 



which translates to the language of A2 (A2 <^ Ai and z ^ 1) as 



Comparing with (^) we get 

/'(I) = ^^^^ (62) 

Now we need to choose the plus sign in all the above formulas as we want a to be positive. 
Otherwise our action will correspond to the solution growing at t ^ —00. That makes us 
choose the first f{z) in (^). 

We would also like to note that according to the equation (|5lD /(—I) = ^/(l) which can 
be checked directly by solving the equations of motion at (f = ir. Moreover, it can be seen 
from (|5l|) that if{—z) is its solution if f{z) is a solution. So we believe there should be 
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some kind of a crossover where the real solution becomes purely imaginary. The fact that the 
correlation function of a real quantity becomes imaginary should not disturb us. It means 
that the correlation function we are computing may not exist for a certain value of Aj and 
can only be understood in the sense of analytic continuation. Apparently, the logarithm of 
probability distribution function which is obtained by a Legendre transform of the action we 
found must remain real. One could perform this transform term by term in our expansion. 
Summarizing everything, the answer for the D-dimensional case is given by 

expjAi }) =exp|yY(pA)5/(cos(/?)| (63) 

while a Legendre transform of the action will give us the probability distribution function in 
the form 



5 iui- 



, 2 u-^ I 9(V33- 1) ,2 1 I 



where ip is the angle between Ui and pi, ■?/' 1, and uj p ^ +oo. 

In conclusion, we would like to say that we have showed by a simple computation that 
WKB calculations are very useful to understand the behavior of randomly driven Burgers 
equation and we hope they will be found useful in other problems of turbulence as well. The 
instanton we found has a spacial homogeneous strain (the velocity was a linear function in the 
inertial interval) and we suspect it to be a general feature of the instantons in the turbulence 
problem. 

We would like to thank A. Polyakov for his valuable advice and useful discussions and V. 
Borue, A. Chekhlov, G. Falkovich, I. Kolokolov, V. Lebedev and V. Yakhot for important 
discussions. One of us (V.G.) is also grateful to J. Maldacena, R. Gopakumar, M. Moriconi 
and C. Nayak for useful conversations. 
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